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Abstract 

An identity is proven that evaluates the determinant of a block tridiagonal matrix 
with (or without) corners as the determinant of the associated transfer matrix (or 
a submatrix of it). 
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1 Introduction 



A tridiagonal matrix with entries given by square matrices is a block tridi- 
agonal matrix; the matrix is banded if off-diagonal blocks are upper or lower 
triangular. Such matrices are of great importance in numerical analysis and 
physics, and to obtain general properties is of great utility. The blocks of 
the inverse matrix of a block tridiagonal matrix can be factored in terms of 
two sets of matrices[10j, and decay rates of their matrix elements have been 
investigated [14J. While the spectral properties of tridiagonal matrices have 
been under study for a long time, those of tridiagonal block matrices are at a 
very initial stage [Tf2] . 



What about determinants? A paper by El-Mikkawyfl] on determinants of 
tridiagonal matrices triggered two interesting generalizations for the evalua- 
tion of determinants of block-tridiagonal and general complex block matrices, 
respectively by Salkuyeh[T5] and Sogabe[T7j. These results encouraged me to 
re-examine a nice identity that I derived in the context of transport [llj, and 
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extend it as a mathematical result for general block-tridiagonal complex ma- 
trices. 

For ordinary tridiagonal matrices, determinants can be evaluated via multi- 
plication of 2 x 2 matrices: 
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Do these procedures generalize to block-tridiagonal matrices? The answer is 
affirmative. If the matrix has corner blocks, the determinant is proportional 
to that of an associated transfer matrix, in general of much smaller size. The 
proof is simple and is given in section 2. A simple modification yields a formula 
for the determinant when corner blocks are absent, and is given in section 3. 
The relation with Salkuyeh's recursion formula is then shown. 



2 The Duality Relation 



Consider the following block-tridiagonal matrix M(z) with blocks Ai, Bi and 
Cj-i (i = that are complex m x m matrices. It is very useful to 

introduce also a complex parameter z in the corner blocks: 
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It is required that off-diagonal blocks are nonsingular: det Bi ^ and det Cj_i 7^ 
for all i. As it will be explained, the matrix is naturally associated with a 
transfer matrix, built as the product of n matrices of size 2m x 2m: 
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where I m is the m x m unit matrix. The transfer matrix is nonsingular, since 

n 

detT= H de1 W lc i-i] ( 5 ) 



The main result, the duality relation, relies on the following lemma: 

(-1)' 



Lemma 1 det M(z) 
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det [ T — z I 2m ] det [ B\ . . . B r , 



Proof: The equation M(z) 1 J r = has a nontrivial solution provided that 
det M(z) = 0, and corresponds to the following linear system in terms of 
the blocks of the matrix and the components ipk G C m of the null vector \1/: 



A 1 xp 1 + B 1 tfj 2 + z- 1 C ^ n = 
Bk^pk+i + A k i) k + Cfc_i^fc-i = 
z B n i\) x + A n ip n + C n _i ^ n -i = 



^ = 2,. ..,n-l) 
The equations (jTJ) are recursive and can be put in the form 
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and iterated. Inclusion of the boundary equations dHj) and jSJ) produces an 
eigenvalue equation for the full transfer matrix (j3j) that involves only the end 
vect or- comp onent s : 











T 


= z 











(9) 



Equation (Q has a nontrivial solution if and only if det[T — zl2 m } — 0, which 
is dual to the condition detM(z) = 0. Both z m det M(z) and det[T — z^m] 
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are polynomials in z of degree 2m and share the same roots, which can- 
not be zero by (J5]). Therefore, the polynomials coincide up to a constant 
of proportionality, which is found by considering the limit case of large z: 
detM(z) « (-l) nm (-z) m det[5i---5„]. ■ 



Before proceeding, let us show that in the special case of tridiagonal matrices 
with corners (m = 1), Lemma 1 with z = 1 yields (pO). 

The factorization 
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is introduced for all factors in the transfer matrix T and produces intermediate 
factors ^ na ^ commute, and allow us to simplify the determinant of the 
lemma: 
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Zi and z 2 are the eigenvalues of the transfer matrix in ([I]) , whose trace is z% + z 2 
and whose determinant is z\z 2 — {b\ - ■ • b n )(co ■ ■ ■ c n _i). ■ 



Multiplication of Lemma 1 by det T 1 gives a variant of it: 
det M(z) = (-l) nm (-z) m det(T- x - -) det[C . . . C„_i] 
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Multiplication of Lemma 1 by the previous equation, with parameter 1/z, 
gives another variant: 



det M(z) det M(l/z) = det 



T + T 



-i 



det[B 1 C ...B n C n . 



Instead of M(z), consider the matrix M(z) — \I nm and the corresponding 
transfer matrix T(A) obtained by replacing the entries Ai with A{ — \I m . Then 
Lemma 1 has a symmetric form, where the roles of eigenvalue and parameter 
exchange between the matrices. For this reason it is called a duality relation. 



Theorem 1 (The Duality Relation) 



det[A/ nm - M{z)} 



-z)- m det[T(A) - zl 2m ] det [Si • • ■ B n ] 



It shows that the parameter z, which enters in M(z) as a boundary term, 
is related to eigenvalues of the matrix T(A) that connects the eigenvector of 
M(z) at the boundaries. 



The duality relation was initially obtained and discussed for Hermitian block 
matrices |ll|ll2lll3j . For n = 2 it is due to Lee and Ioannopoulos[9j. Here I have 
shown that it holds for generic block-tridiagonal matrices, and the proof given 
is even simpler. The introduction of corner values z and 1/z in Hermitian 
tridiagonal matrices (c& = b* k ) was proposed by Hatano and Nelson [7] in a 
model for vortex depinning in superconductors, as a tool to link the decay of 
eigenvectors to the permanence of corresponding eigenvalues on the real axis. 
It has been a subject of intensive research [ToTlolloTnS] . The generalization to 
block matrices is interesting for the study of transport in discrete structures 
such as nanotubes or molecules 18. 3. 191. 



3 Block tridiagonal matrix with no corners 



By a modification of the proof of the lemma, one obtains an identity for the 
determinant of block-tridiagonal matrices M^ -* with no corners (S n = Co = 
in the matrix (j3J)): 

Theorem 2 det M<°> = (-l) nm det[T??] det [Si • ■ • B n _i] 
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where is the upper left block of size m x m of the transfer matrix 
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Proof: The linear system M(°)\I/ = can be translated into the following 
equation, via the transfer matrix technique: 
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Right multiplication by the nonsingular matrix 
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and rewriting the right-hand vector as the product 
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transform (TTTT) into an equation for the transfer matrix T^, that connects the 
boundary components with ^Wi = an d "^o = 0: 






= T (0) 












(12) 



Equation (TTJ) implies that detTjf = 0, which is dual to detM^ = 0. The 
implication translates into an identity by introducing the parameter A and 
comparing the polynomials det[A/ nm — M^ )] and detT( )(A) (obtained by re- 
placing blocks Ai with Ai — XI m ). Since both are polynomials in A of degree 
nm and with the same roots, they must be proportional. Their behaviour for 
large A fixes the constant. ■ 
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For tridiagonal matrices (m = 1) blocks are just scalars and, by means of (fTOj) . 
one shows Theorem 2 simplifies to (f2"1) . 



The formula for the evaluation of det M^ ^ requires n — 1 inversions B^ 1 , mul- 
tiplication of n matrices of size 2m x 2m, and the final evaluation of a deter- 
minant. Salkuyeh[T5] proposed a different procedure for the evaluation of the 
same determinant: 

n 

det M (0) = J] det A fe 

k=l 

Afc = A k — C k _ik k \B k _i, A 1 = A 1 

It requires n — 1 inversions of matrices of size m x m, and the evaluation of 
their determinants. I show that the two procedures are related. 

The transfer matrix T^ ) = T(n) is the product of n matrices. Let T{k) be the 
partial product of k matrices. Then: 



T(k) 



Im 



T(k-1) 



This produces a two-term recurrence relation for blocks 
T(A;)ii = -B^A k T(k - l)n - B k l C k -{T(k - 2) n 



with T(l)n = —B l A\ and T(0)n = I m . The equations by Salkuyeh result 
for A k = -B k T(k) 11 [T(k-l) 11 ] 
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